Retrieving qubit information despite decoherence 
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The time evolution of a qubit, consisting of two single-level quantum dots, is studied in the 
presence of telegraph noise. The dots are connected by two tunneling paths, with an Aharonov- 
Bohm flux enclosed between them. Under special symmetry conditions, which can be achieved by 
tuning gate voltages, there develops partial decoherence: at long times, the off-diagonal element 
of the reduced density matrix (in the basis of the two dot states) approaches a non-zero value, 
generating a circulating current around the loop. The flux dependence of this current contains full 
information on the initial quantum state of the qubit, even at infinite time. Small deviations from 
this symmetry yield a very slow exponential decay towards the fully-decoherent limit. However, 
the amplitudes of this decay also contain the full information on the initial qubit state, measurable 
either via the current or via the occupations of the qubit dots. 

PACS numbers: 03.67.-a, 05.40.-a, 03.65.Yz, 03.67.Pp 



I. INTRODUCTION AND SUMMARY 

Quantum computation operates on information stored 
in "qubits" , which are superpositions of two basic quan- 
tum statesji 

\ip ) = coscv|l) +e l7 sina|2) , (1) 

with two real parameters a and 7. In one realization of 
a solid-state qubit, the two basic states represent single- 
level quantum dots^ where the superposition state (QJ, 
representing a single electron on the two dots, may be 
given as input, or modified by tuning the dot energies 
e 1 2 and the inter-dot tunneling J 12 . In a tight-binding 
language, the Hamiltonian of the qubit is given by 

•H q = e 1 a\a 1 + e 2 a\a 2 - (J 12 a\a 2 + h.c.) , (2) 

where creates an electron on dot n, aJjO) = \n). The 
average energy e = (e 1 + e 2 )/2 does not affect the dy- 
namics of the qubit ^ and therefore we set it equal to 
zero. The dynamics is then determined by the energy 
gap A = e 1 - e 2 and by J 12 , 

n q = (A/2)(a t 1 a 1 - a\a 2 ) - (J 12 a\a 2 + h.c.) . (3) 

In the literature on NMRp- H is often written in the 
equivalent pseudo-spin form 

H q = B • a , (4) 

where <x represents the three Pauli matrices and where 
B z = A/2, B x — iB y = —J\2- Another physical real- 
ization of a qubit involves a superconducting Joscphson 
junction^ 



Clearly, quantum computation requires the stability of 
the quantum state stored on each qubit, and therefore 
it can be used only while this state remains coherent £ 
Interactions between qubits and their environment, in- 
cluding input-output measurement devices, can cause de- 
coherence which destroys the information stored in the 
qubits. Therefore, it is important to study the time 
evolution of the qubit state in the presence of the en- 
vironment. To concentrate on the state of the qubit, 
one traces over the states of the environment, ending up 
with the 2x2 reduced density matrix of the qubit it- 
self, p(t) = Tr cnv [|^(i))(*(<)|], where |*) is the combined 
state of the qubit and the environment. In many cases, 
the coupling to the environment yields full asymptotic 
decoherence, for which the elements of reduced density 
matrix approach the fully-mixed state, 

Pn m (t ^ 00) = S nm /2 , (5) 

independent of which basis is used for the Hilbert space. 
In these cases, the information on the initial qubit quan- 
tum state is totally lost. However, in some symmetric 
cases there exist decoherence-free subspaces, which de- 
couple from the environment, so that at least some of 
the information on the initial quantum state remains 
protected 

Here we discuss special cases in which the full infor- 
mation on the initial qubit state can be retrieved, even 
after a long time, despite decoherence. Since our results 
rely mainly on symmetry, we expect them to hold when- 
ever the required symmetry holds, irrespective of the spe- 
cific nature of the environment and its coupling to the 
qubit. To demonstrate our point we consider the simplest 
model for decoherence, where the environment generates 
a single parameter which fluctuates randomly between 
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two values, as in telegraph noised We thus replace the 
Hamiltonian of the qubit by H q —> H q + f(t)V, where 
f{t) jumps stochastically between +1 and — 1. Indeed, 
such jumps in f(t) may arise e.g. due to equilibrium or 
non-equilibrium sources of noise, e.g. from background 
(natural) charge fluctuations 1 ^*^ or to a capacitive cou- 
pling to a current which flows through a (tunable) neigh- 
boring single-electron transistor] 14 ! 15 

In principle, the Hamiltonian of the coupling between 
the qubit and the noise source, V, may involve the same 
operators which appear in Eq. ([3]): 

V = (CA/2)(4ai - 4« 2 ) - {CAa 2 + h.c.) , (6) 

where the £'s are fixed coefficients, whose size measures 
the coupling between the qubit and the environment. In 
the special cases which we discuss, V commutes with % . 
This requires specific ratios between the coefficients in 
Eq. © and those in Eq. Such ratios can be achieved 
experimentally by tuning gate voltages which control the 
coefficients in Eq. ([3]), or by a careful placing of the 
source of the noise relative to the qubit. When these 
conditions are obeyed, one can switch to a basis which 
diagonalizes H . In this basis, the diagonal elements of 
the reduced density matrix are independent of time, while 
the off-diagonal elements decay to zero, reflecting pure 
dephasing. Translated to the dot basis of the Hilbert 
space, this implies that 

Tr cnv [(*(t)|?g*(i))] 

= y [Pn(t) ~ PzS)] - 2Rc[J 12 p 21 (*)] = co nst. , (7) 

independent of time. Thus, the density matrix never 
reaches the fully- mixed state ([5]). We refer to this sit- 
uation as partial decoherence. 

The decoherence of qubits due to telegraph noise has 
been treated in several earlier papers *2r— In a situation 
where both A and J 12 are present, and both are noisy, 
the qubit 's reduced density matrix usually decays expo- 
nentially towards the fully-mixed stated Alternatively, 
Itakura and Tokura^ considered the special case without 
a gap between the dot energies, A = ( A = [cf. Eqs. 
([21) and JBJ], and found that when J 12 and £j are real, 
then both % and V are symmetric under the interchange 
1 O 2, and therefore the 'bonding' and 'anti-bonding' 
symmetric and antisymmetric states |±) = (|f)±|2))/v / 2 
are eigenstates of both. In that case, the off-diagonal 
element p, oscillates and decays to zero with the 'de- 
phasing time' T 2 . Below we show that this example is 
a special case of a broad family of systems, all of which 
exhibit partial decoherence [cf. Eq. 0]. 

The procedures proposed below require that J 12 should 
be complex and tunable experimentally. To achieve this, 
we connect the two qubit quantum dots via two separate 
tunneling channels, with energies J u and J d (Fig. [I}, 
while a magnetic flux $ is enclosed between them. Uti- 
lizing gauge invariancc, the combined tunneling coupling 



J 12 becomes 

Ju = Ju + J d ^^\Ji2\e 10 , (8) 

with the Aharonov-Bohm phase <j> = 27r$/$ , where 
3?o = hc/e is the flux unit. Both J u and J d (which are 
chosen real) can be tuned via gate voltages, and the phase 
6 can be tuned via the magnetic flux. 




FIG. 1: The two-dot qubit with tunneling channels u and d. 

To demonstrate the result |7j) we present what we term 
'the symmetric case', in which there is no energy gap, 
A = ( A = 0, but J 12 can still be complex. The vanish- 
ing of C A is achieved when the source of noise is located 
symmetrically relative to the two dots, or when the cor- 
relation length of the noise is larger than the qubit's size. 
For an arbitrary flux $, and for given values of the noise 
coefficients ( u and ( d [defined via J u d -> J u d + f(t)( ud ], 
we show that an appropriate tuning of cither J u or J d 
(via corresponding gate voltages) suffices to bring the 
system into a special symmetric case, in which its re- 
duced density matrix never reaches the fully-mixed state 
([5]). Instead, at t — > oo it approaches the asymptotic 
limit 

V2 e i0 Re{e~ ie p 12 (0)] \ 

P ~* \ e- i0 Re[e- ie p 12 (O)} 1/2 J' W 

where p 12 (0) = (I l^oXV'ol 2 ) = e"^ sin(2a)/2 [see Eq. 
dU)] and where 8 is the phase of the complex J 12 , see Eq. 
(|5J). Equation © is a special case of Eq. ([7]), for A = 0. 

The non-zero complex off-diagonal element p 12 gener- 
ates a circulating current around the loop, which can be 
used to retrieve the qubit's information. The equation of 
motion for the number operator a\a 1 is 

d t aW=i u +i d , (io) 

where the operator T u and T d represent the currents into 
dot 1 via the channels u and d, respectively: 

K = i[J»ai<h - h - c -] . 2rf = Vd^A a i ~ h -c-] ■ (11) 
Therefore, the net current into site I is 

X = m± u -X d \y) = (*|2±„ - d t a\am ■ (12) 
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In the asymptotic stationary limit we can drop the last 
term. Tracing over the environment then yields the con- 
ditional average (^(t)\ J u a\a 2 \^> (t)) — > J u p 21 - For the 
symmetric case, we find below that averaging over the 
noise yields 

I — >• 2 J Im[p 12 (oo)] = J sin(2a) sin 9 cos(9 + 7) , (13) 

where J is the average of J u + ( u f(t) over the noise. As 
expected, this current vanishes when <fr is an integer mul- 
tiple of 7r [when also 6 = 0, see Eq. (|SJ]. However, at non- 
trivial fluxes I is non-zero, despite decoherence. This 
current generates an orbital magnetic moment of the elec- 
tron circulating the loop. Measuring its <p— dependence 
can yield both a and 7 [see Eq. ([1}], namely the full in- 
formation stored initially! Unlike the usual equilibrium 
persistent current, which is an odd function of the flux 
(as required by time- reversal symmetry)^ the current 
here is neither odd nor even in the flux. This peculiar 
flux dependence apparently results from the averaging 
over the noise, which breaks time-reversal symmetry. 

Following our analysis of the symmetric case, we also 
consider small deviations from this symmetry, and find 
that such deviations lead asymptotically to the fully- 
mixed state, as also found e. g. in Refs. |12H14| and in ref- 
erences given there. However, for small deviations from 
symmetry there is a distinct separation of time scales. 
After a transient oscillatory stage, the elements of the 
density matrix develop a very slow simple exponential 
decay towards the asymptotic fully-mixed state. These 
slowly decaying terms also include a non-zero difference 
in the dot occupations, z = p xl — p 22 (which approached 
zero relatively quickly in the symmetric case). This dif- 
ference can also be measured experimentally. We show 
that the amplitudes of these exponential terms also con- 
tain the full information on the initial state of the qubit. 
This information can therefore be extracted even after 
the transient stage. 

The plan of the paper is as follows. The formalism 
for the telegraph noise is reviewed in Sec. II. Section III 
then presents several physical environments which can 
generate telegraph noise of the kind discussed here. The 
general conditions for partial decoherence, and the ex- 
ample of the symmetric case, are presented in Sec. IV. 
In Sec. V we then introduce deviations from symmetry, 
and Sec. VI contains a discussion of our results. 



II. TELEGRAPH NOISE 

A treatment of the equation of motion with the under- 
lying stochasticity in f(t) can be found in the literature 
on the theory of Iineshapes4ii2 Here we follow Blume^ 
and average the density matrix p(t) over the histories of 
the stochastic noise, under the condition that at time t 
the random function f(t) has the values b = 1 or — 1. We 
then define a 2— component vector (denoted by bold let- 
ters) pit), such that its b— th component represents this 



conditional average p(t,b). At the end one may average 
over the stochastic process, 

p{t)= £ p{t,b) . (14) 

6=±1 

The function f(t) follows a Markov process^ it jumps 
randomly from 1 to —1 (or from —1 to 1) with the rate 
w , (or w_| ). These jumps in f(t) result from a con- 
tact with some noise source. The noise distribution is 
characterized by the probabilities p ± to find f(t) at the 
values ±1. Detailed balance then implies the relation 

p_w_i = p + w |_, and therefore the jump rates can be 

written as 

W ±T = X P± ' ( 15 ) 

where A = + w ^ represents the inverse time asso- 
ciated with the noise. 

Our main purpose here is to calculate the time evolu- 
tion of p. The equations of motion for the conditional 
averages p(t, b) are 

d tP (t,b) = -i[H q + bV, P (t,b)] 

+ w b _ b p{t, -b) - w_ b;b p(t, b) , (16) 

where we use H = 1 throughout. The first term on the 
right-hand side applies if f(t) remains unchanged at time 
t (i.e. stays equal to b). In this case, the time evolution 
of the density matrix proceeds with the Liouville opera- 
tor which corresponds to the original Hamiltonian, with 
f(t) = b. The last two terms arise if f(t) flips exactly at 
time t, either from —6 to b (second term) or from b to —6 
(last term). 

Each element of the 2x2 reduced density matrix now 
becomes a 2— component vector, p nm , and Eq. (|16|) can 
be written in matrix form, 

*(I^t — W)p m „ = A nm p nm — '^^(J n (Pe m ~ Pnl^lrn) 7 

t 

(17) 

with n,m = 1,2. Here, each parameter in the Hamilto- 
nian Tiq + bV is replaced by a diagonal 2x2 matrix. For 
our specific two-dot system, A m „ = 3 nm = for n = m 
while A 12 = A = AI + ( A cr z represents the energy gap 
variable and J 12 — > J 12 = "A 2 I + C,j<j z represents the 
hopping matrix element (I is the 2x2 unit matrix) . The 
relaxation matrix W consists of the stochastic hopping 
probabilities of the noise, 

W=( ~ W -+ W +- )=X(T- 1) (18) 
V w |_ — iu_| J y ' v ' 

[see Eq. (Ti"5j) ]. where 

T = ( P p + _ P p + _ ) = [I + * x + A P K + ia y )}/2 (19) 
and A.p = p + — p_ . 
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The matrix W determines the time evolution of the 
matrix P(t), where (6|P(t)|a) is the probability for the 
stochastic variable to start at t = with the value a and 
end at t > with the value b: P obeys the equation 
d t P(t) = WP(i), and therefore 



P(t)=e wt = T+(I-T)e- At , 



(20) 



where we have used the identity T 2 = T. At infinite 
time, P approaches T, and thus (&|P(oo)|a) = p b (a, b = 
±), independent of the initial value a. Below we use 
P(£) for some of the solutions for the density matrix. 
It is customary to characterize the noise by its spectral 
function, defined via 

/>oo 

S(u) = 2Re / ^2 Pb , (b-b)(b\P(t)\b')(b' -b)e lU)t dt , 
Jo b,b> 

(21) 



where b = J2bPb b = AP- Using Eq. pO)) one finds 



S(w) 



iw_i w_ + _ 8Xp + p_ 



\{lu 2 + A 2 ) ~ cu 2 + A 2 



(22) 



Below we relate various decay times with special values 
of S(u). 

A convenient way to solve Eqs. (|17[) is by employing 
the Laplace transform, 



poo 

p(s) = / dte- st p(t) . 
Jo 



(23) 



The equations of motion (|17[) then become 

(Is - W)p„ m = p nm (Q) - iA nm p nm 

+i J2( J nePem - PnI J e,n) ■ ( 24 ) 

t 

Since the initial values of the density matrix, p„ m (0) = 
(n|V'o)(V'ol 771 ): do not depend on the stochastic noise, we 
assume the latter to be in its steady state, and identify 
the initial vector p nm (0) with 



,(0) = /w(0)Po, P 



P- 



(25) 



Note that p is an eigenstate of T, with eigenvalue 1, 
and therefore P(£)p = p , so that p corresponds to 
the steady state of the stochastic noise. 

One major issue in this paper concerns the asymptotic 
limit of the density matrix, at long times. Without noise, 
the eigenvalues of H q are ±fi/2, with eigenstates |±). In 
the basis of these eigenstates, p±± remain constant in 
time, while p + _{t) = e -lSlt p+- (0), with the Rabi fre- 
quency 



A 2 



■4|J 12 | 2 



(26) 



In the presence of noise, the density matrix often ap- 
proaches a stationary state, so that d t p — 0. This is 



indeed the case in our analysis. One can then find the 
stationary state by solving the homogeneous linear set of 
equations (| 1 T[) in p nm . Alternatively, one can use Eqs. 
(f2"4"]). together with the identity 



III. 



lim p(t) = lim sp(s) 

t— >oo s— »0 



PHYSICAL MODELS YIELDING 
TELEGRAPH NOISE 



(27) 



The reduced density matrix is obtained by solving the 
equations of motion for the joint density matrix of the 
qubit and the environment, and then tracing over the 
environment degrees of freedom. This procedure is quite 
complicated when the time-dependence of the environ- 
ment degrees of freedom is influenced by those of the 
qubit. This influence is called 'back action'. Neglecting 
this back action implies that one can calculate the time 
dependence of the environment degrees of freedom sepa- 
rately, independent of the qubit states. In the simplest 
model discussed here, the qubit couples to the environ- 
ment only via one degree of freedom, which is represented 
by its time dependent value f(t). Here we review three 
examples of models which have been treated in the liter- 
ature. 

Itakura and Tokurai 2 - already reviewed the literature 
on background charge fluctuations (see also Refs. [ID 
and [l6h . In that case, a single impurity near the qubit 
is either occupied by an electron or empty, with proba- 
bilities p + and p_ and with hopping rates given by Eq. 
(TTSl) . in which A is proportional to e~^ E l k B T anc l AE 1 is 
the activation energy of the impurity (we asume a large 
Coulomb blockade, preventing double occupancy). Ne- 
glecting the back action of the qubit onto the impurity, 
this model reduces to the classical telegraph noise one, 
in which the impurity-qubit coupling generates different 
coefficients in the intra-qubit Hamiltonian for each state 
of the impurity. 

Another possible model concerns a two-level system 
with an energy gap AE (e.g. a double potential well 
created by two neighboring impurities)^ which can be 
represented by a pseudo-spin 1/2. At equilibrium with 
a heat bath at temperature T, the occupation proba- 
bilities obey a Boltzmann distribution p, = 1 — p_ = 
[1 + e~ AE / kB ' T }~ 1 . Each state of this pseudo-spin gener- 
ates different values for the coupling parameters within 
the qubit Hamiltonian, again yielding the telegraph noise 
picture. 

It is usually not easy to justify the neglect of the back 
actions in the above two examples. However, as argued 
by Galperin et al.^ back action may be ignored when 
the dynamics of the fluctuating background charge or the 
two-level system is governed by its coupling to a thermal- 
izing heat bath, which is much stronger than its coupling 
to the qubit. The telegraph noise model is also justi- 
fied in the limit of a very high temperature of this heat 
bathJ£ 
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Here we concentrate on yet another example, in which 
a current between a left and a right reservoirs L and 
R respectively (held at chemical potentials p L > p, R ) 
flows through a single electron transistor (SET), located 
near the qubit. Unlike the above two examples, here the 
fluctuator is not at equilibrium. The states of the en- 
vironment (SET plus reservoirs) include states in which 
an arbitrary number of electrons have moved between 
the two reservoirs, while the SET can be (singly) occu- 
pied or empty. The quantum equations of motion for 
the density matrix of the combined qubit-SET system 
were analyzed in Refs. [ID and [T^. When the bias volt- 
age \i L — p, R is much larger than any other energy in 
the problem, these authors traced over the environment 
states, and obtained equations of motion for the reduced 
qubit density matrix, which are equivalent to our Eqs. 
([16]) . In these equations, the rate of an electron entering 
the SET from the left reservoir, T L , was identified with 

w, , and the rate of an electron leaving the SET to the 

right, r^, was identified with w ,. Here, T L (T R is the 

partial width of the SET level, caused by its coupling to 
the L (R) reservoir. 

When the SET is placed near the qubit, the electron 
on the latter feels an additional Coulomb potential gen- 
erated whenever the SET dot is occupied. Denoting the 
creation operator of an electron on the SET by c , the 
coupling between the SET and the qubit is given by 

?Ant = 4> c o{ U i a l a i + U 2 a\a 2 - [Uja\a 2 + h.c.]) . (28) 

The energy Uj represents a sum of two matrix elements, 
associated with the effect of the SET on the hopping 
between the qubit dots. Assuming the geometry of Fig. 
1, and using the same gauge choice as in Eq. ((5]), these 
matrix elements can be written as 

Uj = U U + U d e» , (29) 

with real U u and U d . Using the conclusions of Refs. [Til 
and [Til, one may replace cjc in Eq. ([28} by a c-number, 
[1 + f(t)]/2. Absorbing the time-independent part in H q 
then yields 

C A = (U t - U 2 )/2 , C u4 = U U J2 . (30) 

These parameters clearly depend on the relative location 
of the SET with respect the two dots and the two tun- 
neling paths. 



IV. PARTIAL DECOHERENCE 

A. General conditions for partial decoherence 

Here we show that a system develops partial decoher- 
ence, i.e. does not approach the fully- mixed state ([5]). 
whenever H and V [Eqs. ^ and ^] commute with 



each other, and discuss the conditions for this to hap- 
pen. The commutator of these operators is given by 

[U v V] = [(Aej-j; 2 ( A )a\a 2 -h.c.} 

+ (Ji 2 0-Ji2Cj)(4a 1 -4a 2 ) . (31) 

This commutator vanishes whenever 

AC./ - JuCa = JuC - JUj = , (32) 

namely 

% = = -3- = K when A 4 , 
A T T* 

^ J 12 J 12 

■Y~ = -yj- = K and C A = when A = , (33) 

J 12 J 12 

where if is a fixed real number. It should be emphasized 
that the conditions (|33|) apply for any environment, when 
f(t) is replaced by an operator acting on the environment 
(e.g. CqCq in the previous section), and are not restricted 
to the telegraph noise example discussed below. 

Equation Q33p requires that Cj/ Ji 2 should be a real 
number. Using the gauge choice ([8]) also for the C's, this 
would require ( U /J u = ( d /J d , or q c = q } , where 

Qj = JjJ d , lQ = UC d - (34) 

Since both J u and J d can be tuned by gate voltages, this 
condition can be achieved experimentally. Alternatively, 
as mentioned in the Introduction, one can imagine con- 
ditions under which £ A = 0, e.g. when the noise source 
is placed symmetrically with respect to the two dots. In 
that case A can be tuned to zero, and Eq. ([32"]) may still 
hold. 

In both of these cases, H q and V can be diagonalized 
simultaneously The corresponding eigenvalues are then 

e ± = ±fi/2, V ±± =±C /2 = ± v /ci/4+|C, / | 2 , (35) 

see Eq. (|26j) . Denoting the corresponding common eigen- 
states by |+) and |— ) (do not confuse with the ± states 
of the fluctuator), the equations of motion in this 
new basis become 

(Is-W)p ±± = p ±± (0) , 

(Is - W)p + _ = p + _ (0) - itlp + , (36) 

with S~2 = S1I + Cfi^z- 

The solution for the diagonal matrix elements is 

p ±± = [Is - WTViitf)) = Pp±±(0) , (37) 

where P is the Laplace transform of P(t), and therefore 

p±±(t) = P(t)p±±(0) = p ±± (0) = p ±± (0) Po , (38) 

where we have used the explicit expressions (|20[) and 
(1231) . The components of these vectors are p±±(t,b) = 
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p b p±±(0), and the factor p b represents the probability 
of finding the stochastic noise at the state b. Averaging 
over this noise [Eq. (fT4|) ]. we find that p ±± (i) = p±±(0). 
Thus, the diagonal elements of the reduced density ma- 
trix do not approach the fully- mixed limit ([5]), implying 
partial decoherence. Note that these matrix elements 
are linear combinations of the original density matrix el- 
ements, with coefficients which involve the mapping from 
1 1,2) to |+,—). In fact, one can repeat the above pro- 
cedure for any operator which commutes with both 'H q 
and V, sec Eq. 0. 

The second equation (j3"B")) yields 



P ±T = F ± p ±T (0) = F ±P±T (0) Po 



(39) 



where 



F ± = [si- W±itl}- 1 

(2s ± 2iQ + A)I + \[a x + Ap{a z + ia y )] T 2i( n cr z 



2[(s ± ifi)(s ± iQ + A) ± i\{ n Ap + <&] 



(40) 



Here we have used Eqs. (|18|) and (|19|) . Averaging over 
the noise [Eq. (|14p]. p H is obtained noting that the av- 
erage of o-^Pq is equal to 1, while the averages of — icr y p a 
and cr z p are equal to Ap. The inverse Laplace transform 
then yields 



P+ _(t)= {A + e a + t + A_e a - t ) P+ _(0) , 



T i 



where 



±(A - 2iC n Ap) + y/A 2 - - 4iAC n Ap 
VA 2 - 4CB - 4iAC n A P 



2a, 



-2if2 - A± 



4iAC n Ap 



(41) 



(42) 



Therefore, p H (t) oscillates and decays to zero. This 

asymptotic limit can also be obtained using Eq. (|2T|) . 
The real parts of — a± represent two decay rates. For 
weak coupling between the qubit and the environment, 
\(\ A, the shorter time (associated with a_) is of or- 
der 1/A, the typical time between the fluctuator jumps. 
To leading order in C/A, the longer relaxation time r is 
approximately given by 



-Re [a , 



4&P+P-/\=&S(0)/2 (43) 



[see Eq. (|22|) ]. This relation with the zero-frequency 
noise spectrum coincides with the well-known 'dephas- 
ing' time T 2 = 2r, generated by fluctuations of the off- 
diagonal coupling between the two energy statesi 4 ' 14 i 23 
Equation (|42j) is also the same as that found (using a dif- 
ferent method) by Itakura and Tokura^i 2 . for the special 
case when A = £ A = Im[J 12 ] = Im[Cj] = 0. 

In the NMR terminology one distinguishes between 
purely dephasing noise, which causes the decay of the 
off-diagonal element of the reduced density matrix (in the 
relevant basis), and purely relaxaional noise, associated 



with the decay of the diagonal matrix elements towards 
1/2. These two types of noise are associated with the de- 
cay times T 2 and T 1 , respectively^^ The model discussed 
in this section describes a purely dephasing mechanism, 
namely T ± — > oo. Similar models have been treated in 
connection with a coupling to a continuum of phoiion 
modest The result (|4"3")l is the same as that obtained for 
white noise, when S(uj) is independent of u. This result 
corresponds to the 'motional narrowing limit' in NMR, 
when the rate of the noise fluctuation A is the shortest 
time in the problem^ 

In the basis of the Hamiltonian eigenstates |±), the 
off-diagonal element p^ decays to zero, and its mag- 
nitude has been used to quantify decoherenceJ^ How- 
ever, this decay is basis-dependent: in the original ba- 
sis of the qubit dot states the off-diagonal matrix ele- 
ments approach non-zero values. This is true whenever 
one encounters partial decoherence. Therefore, using off- 
diagonal elements to characterize decoherence may be 
misleading. It is better to quantify decoherence via a 
basis-independent measure, e.g. [1 — Trp 2 ] = 2(p 11 p 22 — 
|p 12 | 2 ). This quantity approaches 1/2 for the fully-mixed 
state, but is larger than 1/2 for partial decoherence. 



B. Symmetric case 

Equation (|3"3"]l refers to two cases: either the Hamil- 
tonian contains an energy gap A, and an associated 
noise parameter £ A , or both of these variables van- 
ish. Since the results are qualitatively the same in both 
cases, we present explicit expressions for the simpler case 
A = ( A — 0. In this case, fi = 2|J 12 |, £ a = 2|Cj|, and 
the two common eigenstates are easily identified as 



\±) = (\l) T e-* e \2 



(44) 



where 8 is defined in Eq. ([8]) (and simultaneously Qj = 
|£j|e ie , since Cj/^12 i g real).— Substituting J 12 from Eq. 



then yields 



tan6> 



COS ( 



(45) 



with q,j = qQ given in Eq. (fMl) . 

We can now use the results from the previous sub- 
section. For that, we need to map the reduced density 
matrix from our original basis {1,2} to the 'bonding 
antibonding' basis ([4"4"]) and back. Substituting the initial 
conditions 

Z(0) = P++ (Q) - p__(0) = 2Rc[e- ie pi 2 (0)] , 

P+-(0) = Ki(0) - P 22 (0)]/2 + zlm[ e - je pi 2 (0)] (46) 

into Eqs. (|38|) and (|4T|) yields p±±(t) and p H (t), and 

the relations 

z(t) = Pll (t) - p 22 (t) = 2Re[p + _(t)} , 

Pl2 (t) = e ie (am[ P+ _(t)} - [p ++ (t) - p__(t)]/2) (47) 
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yield the reduced density matrix in the original basis. 
It is now easy to check that at long times the diagonal 
elements Pll and p 22 approach 1/2, but the off-diagonal 
element approaches Pi 2 {t —t oo) — > e l6 ~Rc{e~ l6 p 12 (0)], as 
in Eq. ©. 

To present our results graphically, we follow conven- 
tional notations^ and write the reduced density matrix 
in the form 

p=(I + r-<T)/2 , (48) 

where the (real) Bloch vector r = (x, y, z) is defined by 

Pn = (l + *)/2, P22 = (l-z)/2, Pl2 = (x-iy)/2 . 

(49) 



The full thick lines in Fig. [5] show the time evolution of 
the average components of the Bloch vector in the sym- 
metric case, for one example of the parameters. Indeed, 
both the real and the imaginary parts of p 12 approach 
finite limits, while z — > 0. These limits, given by Eq. ©, 
are shown by thin lines. 




t 



FIG. 2: The averages of x, y and z [Eq. (|49])]for cf> = ~.3n, J d = .5, ( d = .3, qj = q c = .5, ( A = 0, p + = p_ = 1/2. The initial 
qubit state [Eq. (JXJ)] is given by a — 7 = .25tt. All energies and inverse times are in units of A. The full (dashed) thick lines 
correspond to A = (0.2). The full (dashed) thin lines represent the exact (approximate) asymptotic behavior. The derivation 
of the dashed lines, for A 7^ 0, is described in Sec. V. 



We now present more details on the derivation of Eq. 
(fl~3|) . The general expression for the circulating cur- 
rent is given in Eq. (fT2"|). With the noise, this equa- 
tion gives the conditional averages, so that we need to 
average over J„p 2 i- I n the- stationary state one has 
Pu -> P12MP0 = e ie Re[e- 8 Vi 2 (0)]Po [see Eq. ®]. 
Substituting also J„ = J u l + ( u er z7 and noting that the 
average of c 2 p is equal to Ap and that J u and C u are 
real, one obtains Eq. (fT5]l . Figure [3] shows the flux- 
dependence of this asymptotic current, for qj = = 1/2 
and several values of the initial qubit relative phase 7 [Eq. 
(PQ)] . Interestingly, the current is odd (even) in <j> for 7 = 
(7 = 7r/2), but is neither odd nor even for intermediate 
values of 7. Equation (|T3)) gives the current in terms of 
the initial qubit parameters a and 7. This current can in 
principle be measured by measuring the orbital magnetic 
moment of the electron on the ring. To retrieve a and 
7, we need to perform three preliminary measurements. 
Fixing the flux at a non-trivial value <p 1 (not an integer 
multiple of ir), one should measure the asymptotic cur- 
rent for two known initial states. These measurements 
determine the device parameters 6 ± and J . Repeating 



the same procedure for another flux (f> 2 and one known 
initial state, one would find 9 2 . The information on an 
unknown initial state can then be extracted by measuring 
the asymptotic current for the same two fluxes. 

V. GENERAL CASE 

Since the calculations for the general case are some- 
what technical, we start with a brief summary of the re- 
sults. Below we present the full solution for the time evo- 
lution of the reduced density matrix, which we have used 
to plot the dashed thick lines in Fig. [5J In this figure, as 
well as in much of the discussion below, we concentrate 
on small deviations from the symmetric case, namely 

I AC,/ - Ji 2 CaI « A, n , \ qj - q c \ < 1 , (50) 

see Eqs. (|32|) - (p4|) . As can be seen from the figure, all 
three components of the Bloch vector [Eq. (|4"9")l ] exhibit 
transient oscillations and then decay with a simple expo- 
nential, 

x Ki x e~ t/T °, y^y e~ t/T °, z w z e~ t/T ° . (51) 
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FIG. 3: The flux dependence of the asymptotic average cur- 
rent (in units of J ) in the symmetric case, for q ; = = 1/2 
and a = .25ir. Increasing dashes correspond to 7 = 0, n/4 
and 7r/2. 
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FIG. 4: The flux dependence of the amplitude z for the same 
parameters as for the dashed lines in Fig. [2l but with 7 as in 
Fig. El 



To leading order in the small parameters ([50|) . the slow 
decay rate is found to be 



4J|(Ci + C^) 



A[A2(fi + C fi Af>)2 + (fi2_ C 2 )2] 



(52) 



where 



n = A + _ = 2Re[J 12 |OI/0] , 
Cn = 2 IOI . -h = (<? ( - 9j) sin ( 

Expanding also in this rate becomes 
, ^ (C n A-C A ^) 2 + 4J|(C|+^) 



(53) 



2ft 2 



S(fi) (54) 



[see Eq. (|22p ] . Again, the relaxation time is related to the 
spectral function of the noise. The corresponding decay 
time, To, indeed becomes infinite in the symmetric limit 
C n A/f2 — ( A = J x = 0, and remains very long for small 
symmetry breaking. This explains the behavior observed 
in Fig. [21 Unlike Eq. (|4~3"]) , the power spectrum function 
S now picks the Rabi frequency f2 of the system, repre- 
senting what Abragam calls the 'adiabatic modulation'^ 
This emphasis on the Rabi frequency is sometimes also 
called the 'rotating wave approximation'. The long relax- 
ation time t can be identified with the relaxational time 
T 1; responsible for the asymptotic decay of the diagonal 
elements of the density matrix towards equal occupations 
(associated with the decay of Z [Eq. (f4"6"f]. which did not 
decay in the symmetric case). 

Below we also evaluate the amplitudes x , y and z 
[Eq. (|5ip]. and the approximate results are shown by 
the thin dashed lines in Fig. [2j The slow exponential 
decay gives an opportunity to measure these amplitudes 
even after a long time. Unlike the symmetric case, where 
asymptotically z — > 0, one now finds non-zero values for 
the occupation difference z = p xl — p 22 - The flux de- 
pendence of the amplitude z is shown in Fig. 2] As in 
Fig. [31 note the even-odd dependence of z (</>). However, 
all these coefficients can be tuned by a few preliminary 



experiments done for a few initial qubit states and a few 
fluxes, as discussed in connection with Fig. [3] 

We now give more details. The unmotivated reader is 
welcome to move to the next section. Since we mainly 
consider small deviations from the symmetric case dis- 
cussed in the previous section, we choose to stay with the 
same basis used for that case, namely Eq. (|44)l . How- 
ever, in order to decrease the number of noise-related 
terms, we choose the phase 9 as the phase of so that 
Cj6 = |Cjl becomes real. This requirement is equiv- 
alent to Eq. (|45|) . upon replacing q } with (in the 
previous section we had q 7 = q^). With the basis (|44p . 
the two Hamiltonian operators become 

H q = (Q/2XI+X+I - I > < I) - (J+-I+X-I +h.c.) , 

V = (Cn/2)(|+><+| - l-X-l) + (Ca/2)(|+X-| + h.c.) , 
J + _ = -A/2 + iJ x (55) 

[see also Eq. ([53")) ]. As a result, Eqs. (|3"o| are now gener- 
alized to the form 

(si - W)Z = Z(0) + 2i(J + _p_ + - J_ + p + _) , 
(si- W±ifl)p ±T = p ±T (0)T*J± T Z . (56) 

Noting that 

Pp = ^(1 + ^0 = ^0, (57) 

the first equation yields 

Z = Z(0)/s + 2iP(3 + _p_ + - J_+p+J . (58) 

Substituting this equation into the equations for p±zp, 
yields two coupled equations for the latter two 2- 
componcnt vectors. 

Solving these equations, and performing the inverse 
Laplace transform, one finds that the time dependence 

of the off-diagonal element of the density matrix (t) 

is a sum over exponential terms, e Q *', where the a.'s are 
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poles of p H (s), found as the roots of a sixth order real 

polynomial, 

6 

d(s) = d t s e . (59) 

1=0 

Substituting the solutions for p±^ into Eq. ([55]) yields 
Z(t). Equations (|47|) are then used to derive the time de- 
pendence of the average Bloch vector (x, y, z), as shown 
by the dashed lines in Fig. [5J Interestingly, there seem to 
be two main time scales. At the beginning one observes 
a transient oscillatory behavior, up to time scales of the 
order given by the symmetric case, Eq. (j42|) . After that, 
all three variables exhibit a very slow pure exponential 
decay, see Eq. (fSTj) . 

To explain this asymptotic behavior, we return to the 
polynomial (|59p . The long-time limit of the solutions 
is determined by the behavior of d(s) at small Laplace 
variable s, where we can use the approximation d(s) « 
<i + d 1 s. The Laplace transform then decays as e a o t = 
e -t / r °, where Tq 1 = — a ~ d /d 1 . This approximation 
is valid as long as d Q /d 1 is small. We find 

d = 4X 2 P+P _ [(C n A - Ca^) 2 + 4 Jf (Ca + Co)] ■ (60) 

Therefore, r becomes infinite, and one has a non-trivial 
stationary solution, only when c?o = 0, which happens 
only when J x = and either A = £ A = or Ca/^ = 
C n /f2 = |C/|/R.c[J 12 e -i6 '], consistent with the results in 

Sec. IV. In all other cases, p H and Z approach finite 

limits as s — > 0, and therefore Eq. (f27|) implies that 
all the components of the Bloch vector approach zero as 
t — > oo, leading to the fully- mixed limit (|5|). 

When the deviation from the symmetric case is small, 
it is appropriate to expand the results in powers of 
(CrjA/fi — £ A ) and J x /Q. In our case, neglecting higher- 
order terms in A, £ A and J x , one has d x « A[A 2 (f2 + 
C a Ap) 2 + (fi 2 - C^) 2 ], yielding Eq. O. 

The coefficients x , y and z in Eq. f|51[) are the 
residues of the poles of the corresponding Laplace trans- 
forms at s — a = — T^ 1 . Therefore, to leading order in 
J ±T , these residues are the same as those for the pole at 
s = in the symmetric case. To this leading order, 

x Q - iy = eT %6 sin(2a) cos(6< + 7) (61) 

and z = [see Eq. ©]. Consequently, the amplitude 
for the slow exponential decay of the circulating current 
is also approximately given by Eq. (|13[) . as plotted in 
Fig. [3] Measuring this amplitude therefore gives the 
same information as discussed in connection with that 
equation (see the end of Sec. IV. B). 

Corrections to the next order in J ±T just shift the val- 
ues of x and y slightly, and therefore we do not discuss 
them here. In contrast, these corrections are crucial for 

z = 2Rc[p_| ] [Eq. (|57])]. since z = at the zeroth 

order. To first order in J H = — (AI + Q A a z )/2 + iJ x l, 

the second Eq. (|56[) becomes 

p + = F + [p + _ (0) - Z(0)/ S ] . (62) 



The amplitude z is given by the average of 
2Rc[lim,^ (sp_| )]. Some algebra then yields 

z « A( Cl C A + c 2 A + c 3 J x )Z{0)/d 1 , 

Cl = X 2 (n + ( n Ap)A P + (O 2 - Ch)(VAp - Co) , 

c 2 = \ 2 (n + ( n A P ) + (n 2 - - ( n A P ) , 

c 3 - 8Ap + p_C n , (63) 

where Z(0) = sin(2a) cos(0+7). The thin dashed lines in 
Fig. [5] were drawn using Eqs. (|51|) . (|6"Tj) and (|6^|) . As can 
be seen, the approximation for z is excellent, while those 
for x and y arc good apart from a small shift which can be 
calculated from the next order. The 0— dependence of z 
is quite complicated, since the coefficients c e also depend 
on 0, via £1 and Qi ■ When qj ^ then J x is proportional 
to sin#, introducing an additional (f>— dependence. 



VI. DISCUSSION 

In this paper we discussed a qubit which is coupled to 
the environment via a single telegraph noise variable fit). 
Apart from the quantum information, which is stored in 
the qubit, the system is characterized by the following 
parameters: the bare energy gap A, the bare hopping 
energies J d and J u = J d qj, and the amplitudes of the 
noise £aj Cd an< 4 C« = Cdlc- To obtain partial decoherence 
we require A = ( A = or C A /A = (/J and q r = q^. 
As already mentioned, the energy gap A and the two 
hopping energies J u and J d can all be tuned by gate 
voltages on the two quantum dots and on the barriers 
along the hopping paths. Therefore, one can in principle 
tune these parameters to the partially decoherent limit. 

Furthermore, £ A depends on the relative locations of 
the noise source and the qubit. To reduce ( A , the noise 
source should be placed symmetrically relative to the two 
dots. In that case, Ci = C2; an d therefore Ca = 0- Also, if 
e 1 and e 2 represent two arbitrary levels of some large dot, 
they will not be strongly affected if the volume of that dot 
is not sensitive to the noise. For an arbitrary qubit state 
and an arbitrary flux, we thus propose to tune A and qj 
until one observes a non-zero asymptotic circulating cur- 
rent, which also generates an orbital magnetic moment. 
After such tuning one can use the same system for re- 
trieving the quantum information, stored initially on the 
qubit, from measuring the current and/or the magnetic 
moment for any other flux and any other initial qubit 
state. 

The measurement of equilibrium persistent currents is 
quite difficult, and it is only recently that novel methods 
were invented to measure themi 2 ^ It remains to be seen 
if such methods can also be applied to the circulating 
currents discussed in the present paper. As mentioned, 
when the system's parameters deviate from the special 
cases with partial decoherence then one can also extract 
the initial qubit information from measurements of the 
occupations of the states on the qubit's quantum dots. 
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Measuring the qubit dots occupations at real time is dif- 
ficult, since the measuring time should be much shorter 
than any decoherence time.— However, in the scheme 
presented here one need not worry about the fast tran- 
sient decay times, and the long decay time r can in fact 
be tuned experimentally. All one needs to do is tune the 
necessary gate voltages and watch for a slow relaxation. 
Therefore, there is a much better chance that existing 
methods for measuring dot occupations will work here. 

Equation ([7]) suggests other options for 'symmetric' 
cases. For example, if J 12 = Cj — then z = p u — p 22 
remains constant in time, while p 12 decays to zero. In 
this case, measuring the time-independent dot occupa- 
tions will yield z(0) = cos(2a) [Eq. JTJ)]. However, to 
gain the flexibility due to the magnetic flux, and to ex- 
tract information on 7, one would still need to deviate 
slightly from this symmetric case. An expansion in J 12 
and/or in Qj would then yield similar slow decays towards 
the fully-mixed state. 



We expect similar qualitative results for more complex 
structures. For example, one can replace each of the 
bonds u and d by a path which goes via a linear chain of 
quantum dots, and one can tune the energy level on one 
or more of these dots through a resonance, thus changing 
the effective hopping energy J u d . As stated, we also 
expect similar results for other sources of noise and for 
a system affected by more than one fluctuator, provided 
one can tune the system to a 'symmetric' limit where the 
commutator [rl q , V] is small. 
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